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Below ~ 630 mK, the ^He atom mass flux, F, that passes through a cell filled with solid hep “^He 
in the pressure range 25.6 - 26.4 bar, rises with falling temperature and at a temperature Td the 
flux drops sharply. The flux above Td has characteristics that are consistent with the presence of 
a bosonic Luttinger liquid. We study F as a function of ®He concentration, x ~ 0.17 — 220 ppm, 
to explore the effect of ^He impurities on the mass flux. We find that the strong reduction of 
the flux is a sharp transition, typically complete within a few mK and a few hundred seconds. 
Modest concentration-dependent hysteresis is present. We hnd that Td is an increasing function of 
X and the Td{x) dependence differs somewhat from the predictions for bulk phase separation for 
Tps vs. X- We conclude that ®He plays an important role in the flux extinction. The dependence 
of F on the solid helium density is also studied. We find that F is sample-dependent, but that 
the temperature dependence of F above Td is universal; data for all samples scales and collapses 
to a universal temperature dependence, independent of ®He concentration or sample history. The 
universal behavior extrapolates to zero flux in the general vicinity of Th ~ 630 mK. With increases 
in temperature, it is possible that a thermally activated process contributes to the degradation of 
the flux. The possibility of the role of disorder and the resulting phase slips as quantum defects on 
one-dimensional conducting pathways is discussed. 

PACS numbers: 67.80.-s, 67.80.B-, 67.80.Bd, 71.10.Pm 


I. INTRODUCTION 

Solid helium is an unique substance that displays a 
combination of classical and quantum properties. It 
has been extensively studied both experimentally and 
theoretically for many decades. One of the most in¬ 
teresting properties of liquid helium is superfluidity, a 
state of matter in three dimensions than occurs below 
a pressure-dependent temperature T\. This quantum 
property is strongly affected by spatial limitation. In the 
two-dimensional (2D) case, the phase transition from the 
superfluid phase to the normal phase is related to the un¬ 
binding of vortices, as described by Berezinskii, Koster- 
litz, and Thouless [i|,0. In one dimension (ID), another 
sort of quantum point defect, the so-called phase slip, 
is responsible for this transition. Quantum Monte-Carlo 
(QMC) simulation Q predicted that the cores of screw 
dislocations in solid helium should be an example of ID 
superfluidity. The flow of superfluid helium in ID can be 
described by the quantum Wdrodynamic theory known 
as Luttinger liquid theory Q. This idea has been con¬ 
firmed by large-scale QMC simulations for the case 
of nanopores. The basic requirements for ID channels 
to demonstrate Luttinger liquid behavior Q are that the 
pore diameter is sufficiently small, the pore length is suf¬ 
ficiently long and the temperature is low enough with 
respect to T\. 

We developed an apparatus, the so-called UMass 
Sandwich^, to study the possible ability of solid helium 
to carry a helium mass flux [Tol - [T^ . Using porous me¬ 
dia, Vycor rods, we are able to apply a chemical potential 
difference between two ends of a solid ^He sample with¬ 
out mechanically squeezing the solid helium lattice itself. 
It was found that an experimental cell filled with solid 


^He can carry a fiuxfl^. but only below some charac¬ 
teristic temperature, Th, and the flux rate substantially 
increases with decreasing temperature. Tiny amounts of 
the impurity ^He also change the flux dramatically [l^ at 
a characteristic low temperature, Td. A brief report that 
discusses some of the ^He concentration dependence has 
appeared [l^. In this report we will describe our measure¬ 
ments as a function of the ^He concentration at several 
pressures in more detail and discuss our interpretations 
of the role of the ^He. 

We note here that this report corrects a thermometry 
error that caused a small shift in the temperature scale 
below ~ 100 mK that was used in the work reported 
previously in Ref. [l^. This was caused by a change 
in the room temperature electronics which we subse¬ 
quently determined introduced a small but measurable 
heating of the thermometer used to measure the solid 
helium temperature. The temperature correction[l^ for 
the work reported in Ref. [Q is less than 1 mk above 
120 mK, 5.5 mK at 80 mK, and can be found from 
TCnew = TC - 0.09637exp(-TC/0.02755). TC is the 
temperature of the thermometer affixed to the experi¬ 
mental cell, Fig. [U and is used to define the solid helium 
temperature, T. All temperatures reported in this work 
include this correction. 


II. EXPERIMENTAL TECHNIQUE 

In this work many freshly grown (and some partially 
annealed) solid ^He-'^He mixture samples have been used 
to study the effect of the ^He impurity concentration, y 
(in the range 0.17 < y < 220 ppm), temperature, and 
pressure on the ^He mass flux. Our experimental meth- 
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ods, have been described in substantial detail in Refs. [ 13 - 
[T3| . We provide a brief discussion of our approach here. 


A. Sample Preparation 


Solid helium samples are grown at constant tempera¬ 
ture from the superfluid in the temperature range 0.3 < 
T < 0.4 K by the condensation of helium into a sam¬ 
ple cell (volume, V = 1.84 cm^) through a direct-access 
heat-sunk capillary followed by an increase in the pres¬ 
sure up to near the melting pressure (about 25.34 bar). 
Subsequent additions of helium are by means of two other 
capillaries in series with Vycor (porous glass with inter¬ 
connected pores of diameter r; 7 nm) rods (0.14 cm dia., 
7.62 cm long). Helium is added to the Vycor to inject 
atoms and create the solid at the desired pressure in the 
range of 25.6 < P < 26.4 bar. A cold plate at the base of 
the sample cell, Fig.[Tl is thermally connected to the mix¬ 
ing chamber of dilution refrigerator and when filled with 
solid '*He can be cooled to about 60 mK. The lowest tem¬ 
perature of the cell is likely limited by the characteristics 
of our 1970’s vintage SHE refrigerator and the heat flux 
through the superfluid-Hlled Vycor rods. Their warmer 
ends have to have much higher temperature than the tem¬ 
perature of the solid-filled cell, up to 1.5 K, to prevent 
the formation of solid helium in the two reservoirs R1 and 
R2 and at the interface between the Vycor rods and the 
reservoirs. The pressure range in the cell has an upper 
limit due to the need to maintain an adequate value of 
the superfluid density in the Vycor so as to not restrict 
the flux of superfluid helium through the Vycor. 

To create samples of known ^He concentration, the cell 
is emptied between each set of measurements [l3|. The 
cell is then filled with nominally pure “^He liquid (0.17 
ppm ^He) up to the saturated vapor pressure through 
use of line 3. Then, a small calibrated volume at room 
temperature is filled with pure ^He to a known pres¬ 
sure. This is injected into the cell via line 3 and this 
injection is followed by additional '*He, which also en¬ 
ters through line 3, to bring the cell close to the melting 
curve. With knowledge of the relevant volumes and pres¬ 
sures, a known concentration of ^He is thus introduced 
into the cell. The solid is then grown by injection of '^He 
through the two Vycor rods. After the sample is grown, 
it is allowed to rest for Ri 5 — 10 hours at a solid helium 
temperature < 0.4 K before starting any measurements. 
Most solid helium samples are freshly grown (and not 
annealed above 0.5 K). As we will see, we find repro¬ 
ducibility in a given sample with temperature changes, 
which suggests that the samples are adequately in equi¬ 
librium after being created. As has been seen previously 
for nominal-purity well helium (measured for this work 
to be ^ 0.17 ppm ^He concentration), high temperature 
annealing leads, on cooling, either to a substantial flux 
decrease or to complete flux extinction, with in that case 
no evidence for flux at lower temperatures. 



FIG. 1. Schematic diagram of the cell used for flow exper¬ 
iments. Two capillaries, 1 and 2, go to liquid reservoirs R1 
and R2 at the top ends of the Vycor rods, VI and V2. Cap¬ 
illary 3 enters from the side of the cell and is used for adding 
helium to the cell. Two capacitance pressure gauges. Cl and 
C2, are located on either end of the cell for in situ pressure 
measurements of the solid '^He. Pressures in the lines 1 and 
2 are read by pressure gauges, PI and P2, outside the cryo¬ 
stat. Each reservoir has a heater, HI, H2, which prevents 
the liquid in it from freezing and allows the temperatures of 
the reservoirs to be controlled. The relevant temperatures are 
read by calibrated carbon resistance thermometers Tl, T2 and 
TC. [Reproduced from figure 1 in Ref. El] 


B. Measurement Procedure 

To initiate the flux, an initial chemical potential differ¬ 
ence, Ap,o, is applied between the tops of the Vycor rods 
by changing the temperatures, Tl and T2, of the two 
reservoirs, R1 and R2, to create a temperature difference 
between them (See Fig.[T]). This creates a flux due to the 
fountain effect, which is seen by observing changes in PI 
and P2. Since we monitor the pressures and the temper¬ 
atures as a function of time, the chemical potential A/i 
can be calculated. 



where p and S are the temperature-dependent density 
and entropy of liquid helium, respectively. In contrast to 
some of the earlier work from our lab, where AP = PI — 
P2 was applied El by direct mass injection, our current 
study uses the application of a temperature difference 
AT = Tl — T 2 E 3 . This approach offers two advantages. 
It allows for smaller density changes in the solid helium 
than was the case for direct injection of “^He to the sample 
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cell through one of the lines, 1 or 2. And, it allows us to 
keep constant the total amount of ^He in the apparatus. 

An example of the procedure used for the flux mea¬ 
surements for a solid "^He sample with a 10.2 ppm ^He 
impurity content is shown in Fig. [51 The creation of a 
change in the energy deposited in heaters HI and H2 re¬ 
sults in a temperature difference, AT = T1—T2, between 
the reservoirs Fig. El^a), R1 and R2, at the tops of the 
Vycor rods and results in pressure responses Fig. El^b), 
PI, P2 and AP = PI — P2 due to the fountain effect at 
a sequence of rising solid helium temperatures Fig. He), 
TC. The derivative of AP, 


F = 


d{Pl - P2) 
dt 


( 2 ) 


is taken to be reasonably proportional to the flux, P, 
of atoms that move from one reservoir to the other. We 
report the rate of pressure change in mbar/s units, where 
0.1 mbar/s corresponds to a flux of Ri 4.8 x 10“® g/s. 

We use measurements of F of this sort in two related 
ways. In one, we study how the flux, P, depends on the 
chemical potential, A/i, as time evolves, as the chemical 
potential changes from its initially imposed peak value, 
A/io, imposed by the initial AT, to zero as equilibrium is 
restored by the creation of a fountain effect induced pres¬ 
sure difference, PI — P2. In the other, we study, for a 
given value of the imposed AT, how the maximum result¬ 
ing flux, P, depends on the temperature of the solid "^He. 
Data, of the sort shown in Fig. [2] for a specific 10.2 ppm 
^He sample, is taken for a variety of solid helium sam¬ 
ples, each with a specific value of the ^He concentration. 
We choose to focus on the behavior of AP, which also 
allows us to eliminate a small long term drift in PI and 
P2, which is typically present in our long-duration mea¬ 
surements due to main helium bath level changes. Our 
basic conclusions are not changed if instead we focus on 
the individual behaviors of PI or P2. 



- 0.01 < 
0.01 


FIG. 2. (color online). An example of flux measurements for 
10.2 ppm sample. Here (a) temperatures are established for 
each of the reservoirs, held constant while the pressure in each 
of the reservoirs stabilizes, and then the temperature values 
of the two reservoirs are interchanged. The interchange re¬ 
sults in a flux of atoms driven by the fountain effect which is 
recorded by (b) the pressure gauges PI and P2. The rate of 
change of PI — P2 provides a measure of the time dependent 
flux, (c) The solid ^He temperature is changed and the pro¬ 
cess continues for a sequence of solid ^He temperatures, TC. 


where A and b are fit parameters. The parameter b is 
temperature independent and as we will seeAn the pres¬ 
sure range of our study b is less than 0.5 [l^, but does 
depend on pressure. We will return to a discussion of the 
characteristics of A and b later. 


D. Flux Dependence on Temperature 


C. Flux Dependence on A/i 

Using data, including that in Fig. [51 F vs. A/t is ob¬ 
tained for positive AT values and presented in Fig. O 
for a set of solid helium temperatures. The maximum 
flux values are typically constrained by the solid helium 
sample. But, for the lower temperatures, the constraint 
is imposed by the temperature of the reservoir at the 
upper end of the Vycor as shown previously!!^. The 
dashed line in Fig. [3] represents the flux limit imposed 
by the Vycor. Similar behavior is seen for negative AT 
values. As was found earlier [l^, a power law provides 
a good two-parameter characterization for data of this 
sort: 

F = A(A/.)^ (3) 


The maximum flux measured through the solid he¬ 
lium that results from a specific imposed AT (typically 
±10mK, Fig. [5]) occurs for a resulting A/t in the range 5-8 
mJ/g, and has a temperature dependence as illustrated 
in Fig. |T| for the case of a different solid sample with 
a ^He concentration of 19.5 ppm. As we will see, this 
general behavior is present for all of the concentrations 
we have studied. This maximum flux, T, increases with 
falling temperature, with warming and cooling showing 
the same values of the flux for a given sample, so long 
as the sample is not annealed. And, there is a sharp re¬ 
versible decrease of the flux at a concentration-dependent 
temperature, T^. 

As seen in Fig. [4] for this x = 19-5 Ppm ^He impu¬ 
rity sample, the flux that results from a given AT is an 
increasing function of decreasing temperature until the 
temperature drops below ~ 105 mK, below which there is 
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FIG. 3. (color online). An example of the F{Afj,) depen¬ 
dence for different solid helium temperatures for the case of 
a sample with a ®He concentration of 10.2 ppm. Solid lines 
are power law fits by use of Eq. The dashed line is the 
upper limit of the flux due to the Vycor bottle-neck for reser¬ 
voir temperatures of 1.48 K (see Fig.3 in Ref. 0)- Note, 
before taking a derivative to calculate F, a moving average 
of the AP{t) data, FigO was determined: by 3 points for 
TC < 0.25 K, by 7 points for 0.25 < TC < 0.40 K, by 9 
points for TC — 0.445 K and by 12 points for TC = 0.492 K. 


no flux. As an illustration of just how sharp and prompt 
the extinction behavior is, consider Fig. [S] and Fig. [51 
These figures illustrate that for a x = 10.2 ppm sample, 
the transition from flux to no flux is no more than ss 1.5 
mK wide, with the cessation of the flux complete within 
no more than ^ 350 seconds. Similarly, in Fig. [5] we see 
that there is no flow at a cell temperature of 99 mK, but 
that an increase in the cell temperature to a fixed value 
near 100 mK results in a growth of the flux, with a flux 
recovery time of ^ 600 seconds. The difference in the 
temperature of the sharp change in F between cooling 
and warming shows a small hysteresis at this value of the 
^He concentration. The sharp gradient in the slope of F 
vs. T near 100 mK seen in Fig.|4]is stable. That is, if the 
temperature remains fixed, then the value of F remains 
stable at any point in the Td transition region. 

Now, the time noted for the flux to make the no-flow 
to flow recovery (or the reverse) likely places some con¬ 
straints on scenarios for what causes the transition from 
a state of no flow to a state of flow. One possibility is 
that it takes this long for the temperature of the solid 
to change. To explore this, a calculation for our cylin¬ 
drical geometry that incorporates the Kapitza resistance 
between the solid and the copper wall and the properties 
of solid helium has been done by Mullinfl^ with the re¬ 
sult that the time for this thermal equilibration to take 
place is predicted to be no more than ~ 50 msec. This 



FIG. 4. (color online). Maximum values for the flux as a 
function of temperature for the case of a solid with 19.5 ppm 
^He at P = 26.40 bar. The sharp behavior of the flux extinc¬ 
tion in a very narrow range of temperature near T — Td is 
evident. 



time (s) 

FIG. 5. (color online). An example of the extinction of the 
flux as the temperature of the solid helium falls below T — Td 
for a 10.2 ppm ®He sample at P = 26.30 bar. This figure 
has the corrected temperature scale and is a revision of the 
similar figure presented in Ref. [l3| 


result is consistent with the thermal e quili bration experi¬ 
ments carried out by Huber and Maris [l^ that indicated 
that equilibration near 100 mK is achieved in ^ 10 msec. 
These facts indicate that the equilibration time for tem¬ 
perature of the solid is much faster than the observed 
recovery times and thus the flux change must be related 
to the movement of the ^He in the solid. We will discuss 
this further later. 


AP (bar) 
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TABLE I. Sample Characteristics 



FIG. 6. (color online). An example of the recovery of the flux 
near T = Td for the same 10.2 ppm sample as in Fig. [H No 
flux is seen in the presence of an imposed AT until the cell 
temperature increases above about 100 mK, after which the 
flux recovers in a few hundred seconds. This figure has the 
corrected temperature scale and is a revision of the similar 
figure presented in Ref. [I^ 


X <5x P (bar) Td (mK) STd (mK) 


0.17 


25.64 

72.5 

7 

0.17 

- 

25.90 

72.5 

3 

1.0 

0.2 

25.86 

80.5 

5 

2.0 

0.2 

26.10 

88.5 

5 

4.0 

0.5 

26.09 

91.5 

5 

10.2 

0.5 

26.30 

97.0 

5 

15.0 

3 

25.92 

99.5 

3 

19.5 

1 

26.40 

103 

3 

25.5 

1.2 

26.12 

106 

3 

40.0* 

5 

26.15 

109 

2 

61.0 

3 

26.36 

111 

2 

119.3 

6 

26.40 

115.5 

2 

220.0* 

30 

25.90 

125 

2 


The ^He concentration (ppm) is in two cases estimated {*) 
based on the ln{x) vs. 1/T linear dependence shown in 
Fig |13l inset. The quantities Sx and STd represent 
uncertainties in the determination of x ^-nd Td- 


E. Flux Dependence on ^He Concentration 

The same procedure shown in Fig. [5] has been used for 
a substantial set of solid helium samples with different 
^He impurity concentrations x that ranged from a low 
for nominal well-helium (again, here measured to be 0.17 
ppm ^He) to a high of 220 ppm, as listed in Table I. Two 
examples of the mass flux temperature dependencies are 
shown in Figs. [7] and |S] for x = 0.17 and 2 ppm, re¬ 
spectively. Data points here represent the maximum flux 
values normalized to the maximum flux rate at 200 mK 
to facilitate the comparison. The relevance of such nor¬ 
malization will become more apparent shortly. 

Figs. [7] and m document hysteresis in the vicinity of Td 
with the flux drop during cooling typically taking place 
at slightly lower temperature than the flux rise during 
warming. This hysteresis is measurable at most of the 
^He concentrations studied and might be considered as 
a feature of a first order phase transition. This suggests 
that phase separation my be important and we will re¬ 
turn to this point later. The hysteresis is most evident at 
the lowest concentrations. The width of this hysteresis 
for low concentrations is shown in Fig. [HI Figs. [7] and [5] 
demonstrate that for low concentrations the flux does not 
drop to zero and recovers on cooling below T < Td- The 
recovery of the flux as the temperature is lowered below 
TwfT^. suggests that the role played by the ^He sat¬ 
urates. We will return to this point in the Comments 
section, section IV. The dashed curves on the figures 
are vertically shifted continuations of the solid smooth 
curves, which serve to characterize the data. The sig- 



FIG. 7. (color online). An example of flux reduction and hys¬ 
teresis near T = Td for a 0.17 ppm ®He and P — 25.64 bar 
sample. The smooth curve will be discussed in the next sec¬ 
tion. 


nificance of these smooth curves will be discussed in the 
next section. If x is more than about 10 ppm, e.g. as 
shown in Fig.|4l then there is no flux recovery down to ~ 
60 mK (the lowest temperature for these measurements). 

Data for a range of samples with different concentra¬ 
tions and sample histories are shown in Fig. [TOl In each 
case the maximum flux value shown is that which results 
from the same value of the imposed AT. The shift in 
Td with concentration is evident. Different samples with 
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FIG. 8. (color online). An example of flux reduction and 
hysteresis near T — Td for a 2.0 ppm ®He and P = 26.10 bar 
sample. The smooth curve will be discussed in the next sec¬ 
tion. 



FIG. 9. (color online). The width of the hysteresis seen in 
the vicinity olT — Td for low ®He concentrations. The width 
of the hysteresis region narrows steeply with increasing con¬ 
centration. 


different histories at a given concentration have some¬ 
what different absolute values of F, but the temperature 
dependence and value of Td are reproducible for a given 
concentration. 


III. DISCUSSION 

A. Universal Temperature Dependence 


We now present in more detail the temperature depen¬ 
dence for temperatures above Td- As we will show, the 



0.17 ppm warming 
0.17 ppm cooling 
0.17 ppm warming 
4 ppm warming 
4 ppm warming 
10 ppm cooling 
10 ppm warming 
19.5 ppm cooling 
19.5 ppm warming 
61 ppm cooling 
61 ppm warming 
119 ppm cooling 
119 ppm warming 
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Temperature (K) 


FIG. 10. (color online). The temperature dependence of the 
flux observed for “^He with several ®He impurity concentra¬ 
tions and experimental conditions, determined in each case 
with a constant value of AT, which yields a maximum flux, 
F, that appears for A/r in the range 5-8 mj/g. This figure 
has the corrected temperature scale and is a revision of the 
similar figure. Fig. 1, presented in Ref. [l^ . For each data 
set the solid pressure was in the range 26 ± 0.4 bar. 


temperature dependence for T > Td is robust, but the 
absolute value of the flux depends importantly on the 
sample and its history. As an example of the sort of vari¬ 
ability that we have found consider the data shown in 
Fig. HU for the sample with y = 19.5 ppm. The flux val¬ 
ues are presented here in mbar/s units (not normalized 
values) in order to compare the behavior of the flux be¬ 
fore and after the temperature was increased to 620 mK, 
where the flux was no longer measurable. When cooled, 
after the sample cell was warmed, the flux was greatly 
reduced (circles). The data set for the larger values of F 
shown here is the same data set shown in Fig. 0] 

What is not immediately apparent in Fig. [10] and 
Fig. [TU is that the temperature dependence at tempera¬ 
tures above the peak flux reached is robust. To demon¬ 
strate this most clearly, we normalize the many data sets 
shown in Fig. (TUI We accomplish this by use of a mul¬ 
tiplicative factor for each data set to force the various 
values of F to superimpose at T = 0.2 K. The normal¬ 
ized flux temperature dependencies for the samples of 
different ^He concentrations are presented in Fig. IT^ 
One can see again here, as was evident in Fig. (TU] that 
the Td values shift to higher temperatures with higher y 
values. At the same time, F{T) for different samples at 
T > Td collapse to a universal temperature dependence 
[l^ . This Figure also shows that in the temperature 
range in the vicinity of the peak value of the flux (near 
T = Td) the peak becomes more rounded with less cur¬ 
vature for larger y values, as does the slope of Td vs. 
X- 
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FIG. 11. (color online). An example of flux extinction near 
T = Td for a 19.5 ppm ^He sample at P = 26.40 bar sample. 
For comparison, data for cooling is shown (solid circles) after 
the flux ceased at T — 0.625 K. Note, flux values are shown 
in mbar/s units. 



0.1 0.2 0.3 0.4 0.5 0.6 

Temperature (K) 


FIG. 12. (color online). The temperature dependence of the 
normalized flux observed for ^He with several ®He impurity 
concentrations and experimental conditions, with the solid 
'^He pressure 26 ± 0.2 bar. Fitted line: see text. 


Td values determined from the data in Fig. [12] and 
other data like them are shown in Fig. [T3| to show how 
the ^He concentration affects Td. It is natural to compare 
these temperature-dependent data to the data on phase 
separation in solid helium. According to Ref. [^, ex¬ 
trapolated to 26 bar, the temperature of Solid-Solid (bcc 
^He-rich inclusions form inside the hep ^He-rich matrix) 
phase separation temperature, Tpg, can be found from 


Tp, = [0.80(1 - 2x) + 0.14] / ln(l/x - 1) (4) 

and this is represented in Fig.[T3|by the dashed line. The 
number 0.80 in this expression comes from the extrapo¬ 
lation to the pressure of our experiment. In the case of 
bulk phase separation for our pressure range, another sit¬ 
uation is present: liquid ^He-rich regions form inside the 
solid '^He matrix (the so-called, Solid-Liquid case). This 
scenario was calculated in Ref.jl^ and is shown by the 
solid line in the same Figure. It can be seen here that 
our Td temperatures lie above the Solid-Liquid case. This 
will be discussed further in Section IIII Cl 



FIG. 13. (color online). Temperature of the sharp drop in 
F, Td- Inset: ln{x) vs. l/T; see text. This figure has the 
corrected temperature scale and is a revision of the similar 
figure, Fig. 5, presented in Ref. [l^ . 

In an attempt to further characterize the data we have 
utilized several functions. In recent presentations [T^ [l3l| . 
we have favored F = A — B ex-p{—E/T), which is moti¬ 
vated by the thought that some thermally activated pro¬ 
cess may be relevant. Independent fits of this functional 
dependence to all of the data sets results in the recogni¬ 
tion of the universal character of the temperature depen¬ 
dence. We find a good characterization of the data with 


F = Fo[l-1.21exp(-F;/r)], (5) 

The value of E that results from such a characteriza¬ 
tion of the data depends weakly on pressure as shown in 
FigHH The higher the pressure (density), the lower the 
value of E, i.e. the F{T) dependence gets steeper with 
pressure. The use of colors for symbols in this Figure is 


































FIG. 14. (color online). Pressure dependence of the parame¬ 
ter E. 

the same as in Fig. [TT] and Fig. [T^l There is no apparent 
dependence of i? on %. 

The value of Fq may be interpreted to be proportional 
to the number of conducting pathways inside the solid 
helium. Le. at T ~ T/j conducting pathways are being 
partially annealed (completely in some cases) leading to 
a substantial flux decrease in the whole Td < T < Th 
temperature range on subsequent cooling. 

As we have pointed out earlier the entire temper¬ 
ature dependence is not fully explained by thermal ac¬ 
tivation since the flux extrapolates to zero at a finite 
temperature. So, whatever controls the decrease in flux 
with increasing temperature must have an explanation 
that goes beyond simple thermal activation. 

As an alternate approach to characterize the tem¬ 
perature dependence, the normalized universal data for 
T > Td from Fig. [12] can be inverted, to obtain 

something we might call a flux resistance as shown in 
Fig. US) This approach allows us to explore whether there 
might be any power law behavior, although the tempera¬ 
ture range is very narrow for such an approach. One can 
see in Fig. [13] that there appears to be a crossover in the 
behavior of the temperature dependence. The data for 
the range of samples studied can be described reasonably 
well by 

(F)-i = F{0.2K)/F = AT^ + BT^, (6) 

where A, B, k and m are parameters. We find that with 
the choice of fc = 1 a fit to the data yields m = 5.8 ± 0.3. 

It is an open question as to what the origin of those two 
apparently distinct contributions to the temperature de¬ 
pendence of the mass flux resistance is. But, such behav¬ 
ior is not without precedent for a quasi-one-dimensional 
system. Consider, for example, the case of supercon¬ 
ducting nanowires below the transition temperature, Tc- 
These nanowires demonstrate nonzero resistance at any 


FIG. 15. (Color online) Temperature dependence of the flux 
resistance, {F)~^, measured through the solid sample (see 
Ref.fl^, Fig.4) and presented here on log-log scale. The solid 
line is a fit of the data by Eq. and the dashed and dotted 
lines represent linear and T®'® behavior, respectively. 

finite temperature, apparently due to the presence of 
phase slips in the order parameter that result in dissi¬ 
pation, which destroys superconductivity. These phase 
slips are due to thermal fluctuations at higher tempera¬ 
tures close to, but below Tc, or to quantum-mechanical 
tunneling at low temperatures, so-called quantum phase 
slips (QPS). Electrical transport measurements in single¬ 
crystal Sn nanowires [2lj and its analysis showed a 
power-law dependence p{T) T“ at T < Tc with an ex¬ 
ponent Of sa 5 for nanowires of diameter 20 and 40 nm, 
but much larger values of a for larger wire diameters. 
These data were interpreted in terms of the unbinding 
of quantu m p hase slips with temperature. Also it was 
predicted [23l| that in the limit of very thin wires and low 
temperatures, where unbound QPS behave as a gas, the 
temperature dependence of the wire resistivity should be¬ 
come linear at the transition to the disordered (i.e. non¬ 
superconducting) phase. As shown in Fig.[T5l the data is 
consistent with a linear dependence of the flux resistance 
for the flux that we observe at low temperature. The rela¬ 
tionship that might exist between these rather different 
physical systems, conducting pathways in solid helium 
and very thin wires, has not been explored theoretically, 
but both systems may be describable by Luttinger liquid 
theory. 


B. Luttinger Liquid 

The non-linear behavior of F vs. Afi shown in Fig. O 
is reasonably well represented by Eq. (3), where b is less 
than 0.5, independent of T. An example of the indepen¬ 
dence of b from temperature is shown in Fig. |TH| for x 
= 10.2 ppm. Here, as in our earlier measurements for 
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nominal purity well helium [l^. all of the temperature 
dependence is contained in the amplitude A. We have 
previously shown that for well helium b depends on pres¬ 
sure. As we have suggested, this non-linear behavior and 
independence of b from temperature supports the possi¬ 
bility that the flux is carried by one-dimensional paths, 
e.g. perhaps the cores of edge dislocations jl^, and can 
be described by the properties of a, so-called. Bosonic 
Luttinger liquid 0. 




0.4 0.6 

SP = P - Pmc( 


FIG. 17. (color online). Pressure dependence of the fit param¬ 
eter b. Open square data points here correspond to nominally 
pure ^He samples (0.17 ppm ®He) and other data points rep¬ 
resent the data for ®He-"^He mixtures with x > 0-17 ppm. 
Here, SP is the distance above the melting curve in bar. 


This suggests that for T > Td there is no measurable ^He 
role in either the F(A^) or the F{T) dependencies. 


FIG. 16. (color online). Temperature dependence of fit pa¬ 
rameters A (a) and b (b) for the data in FigO see Eq.Q. 


As noted, the exponent b is temperature independent, 
but it depends on the solid helium pressure. The higher 
the pressure, the larger is the value of b. Data for b as 
a function of the distance from the melting curve, 5P = 
P — Pmc is presented in Fig. [T7] for nominally pure (170 
ppb) helium as well as for a number of concentrations. 
Although the number of concentrations for which we have 
data adequate to determine b for a range of pressures is 
limited, there is apparently no significant dependence of 
b on X. 

The data for b above can be used to obtain the pres¬ 
sure dependence of the Luttinger parameter, g. If we 
presume that we have a number of independent random 
scattering sites that introduce phase slips, then the Lut¬ 
tinger parameter ocan be obtained from b by means of 
g = [(1/6) -I- l]/2 1^. The results shown in Fig. fTSl sug- 
gest that for such a scenario we are well in the Luttinger 
regime, but that with increasing pressure g decreases and 
we approach the non-superfluid regime. This is consis¬ 
tent with previous work[l^ [lll| in which the flux dis¬ 
appeared at higher pressures. The colors of symbols in 
Fig. HH] have the same sense as in Fig. [TT] 

Based on the pressure (density) dependence of g and 
E, FigsdH] and [TTl respectively, one can see that ^He 
impurities in its range studied do not affect these data. 
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FIG. 18. (color online). Pressure dependence of the Luttinger 
parameter, g, presuming that g = [(1/6) -I- l)/2. Again here, 
SP is the distance above the melting curve in bar. 


C. COMMENTS 

We have discussed the effect of ^He impurities on the 
flux measured and have shown that Tpg ^ Td- We sum¬ 
marize the ^He effects that we have observed. First of 
all, there is a sharp flux extinction at a characteristic 
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temperature Td, which itself depends on ^He concentra¬ 
tion, X, and this Td{x) dependence (see FiglT^ is rem¬ 
iniscent of the bulk phase separation shifted to higher 
temperatures. Another point is that this extinction is 
a slow process compared with the thermal equilibration 
time for the solid lO’s of msec, but a fast process 
(~ lOO’s of seconds), compared to the time required for 
complete solid phase separation (dozens of hours HI) 
in samples of higher concentrations than we have used 
here. These effects could suggest that only a small frac¬ 
tion of ^He is responsible for the flux extinction and this 
phenomenon is due to ^He redistribution in the vicinity 
of the phase separation temperature. Work by Edwards 
et al.[^ indicates deviations from behavior in the 
specific heat for T > Tp.,, which suggests local ^He con¬ 
centration fluctuations [2 7j| are likely relevant to our ob¬ 
servations; the ^He fluctuates in position and is thereby 
able to block the flux at a temperature above the bulk 
phase separation temperature. It is the case that for the 
very low concentrations we have for the most part used 
in this work, diffusion can be quite fast. We will explore 
this further below. 

Based on the F ~ (A/x)^ dependence, as was already 
shown in Ref. [T^ . we suggested that the flux could be 
consistent with a one-dimensional scenario. One can en¬ 
vision two possible candidates for these ID pathways: (1) 
liquid channels, e.g. between grain boundaries and the 
sample cell wall[^ or (2) the cores of dislocations in solid 
helium. 


1. Liquid Channels 

Although we have previously argued that liquid chan¬ 
nels are likely not present, we discuss them further here. 
Were they present, the superfluid density in them could 
be reduced or eliminated by the migration of ^He. Un¬ 
fortunately, there is no reasonable estimate of how many 
of these liquid channels there might be. In the pressure 
range we have studied, the diameter of these channels 
can be calculated to be 6 — 31 nm; the higher the 
pressure, the smaller the channel diameter. In the most 
extreme example, if one such liquid channel were to span 
the distance between the Vyror rods with diameter of 20 
nm, to fill it with ^He would require 1.4 x 10^^ atoms. At 
10 ppm ^He concentration there are enough ^He atoms 
present in our solid mixture samaple to fill ~ 3 x 10® 
such channels. Although we don’t have a good estimate 
of the number of such channels, given the expected dif¬ 
fusion times it appears that diffusion to such channels 
would be a fast process. 

In work reported in Ref. [2^, torsional oscillator, 
TO, measurements to determine the superfluid density, 
Ps, were carried out for nanometer-size channels (folded 
sheet mesoporous materials) of diameter D = 1.5 — 
4.7 nm (^. Isojl filled with superfluid helium. This study 
revealed a transition from a Kosterlitz-Thouless behavior 
to a ID-like temperature dependence of the apparent su¬ 


perfluid density only for D < 2.2 nm. The temperature 
dependence they found, ps for D — 1.8 nm [2^, can be 
well fit by Eq.® which has been chosen to fit our flux 
temperature dependence, F{T) (see Fig. [Tn]). A value of 
E ^ 0.4 K for the parameter E is found for the data 
of Ref.H^. This functional dependence, which is not 
present in the work of Ref. [2^ above D = 2 nm, lends 
support to the notion of a ID scenario for our observed 
flux and further suggests that the liquid channels that are 
predicted would be too large in diameter to demonstrate 
1-D behavior. 



T(K) 


FIG. 19. (color online). Temperature dependence of the tor¬ 
sional oscillator data from Ref.[2^ for a pore diameter of 1.8 
nm. Here Af/Af{T = 0) is the relative frequency shift seen 
in the torsional oscillator as a function of temperature. The 
smooth curve is a fit of the data directly to Eq. that we 
have used to characterize the universal temperature depen¬ 
dence of our flux through solid helium; we find in this case 
E ~ 0.4 K. If we convert to a flux resistance in this case, the 
power law, Eq. ®, does not provide a reasonable fit. 


2. Dislocations 

Based on the possibility of ID of conducting 
pathways [HI, one approach is to assume that ®He im¬ 
purities bind on dislocation cores (or their intersections) 
in the solid helium and block the flux that is carried by 
such pathways. This notion is supported by QMC sim¬ 
ulations [ 3 ^ that show that ®He impurities diminish the 
superfluid density along the core of screw dislocations in 
hep ^He by binding on them. To illustrate this point, the 
inset to FiglTS] shows ln(x) vs 1/T. Straight lines here 
are presented by y = exp(—i?/T) with R approximately 
independent of temperature, and R = 0.94 K and 1.02 K 
for Solid-Solid (dashed curve) and Solid-Liquid (solid 
curve) bulk phase separation, respectivelyfl^. A fit of 
the Td data (squares, dark solid line) by y = exp(—i?/T) 
gives R = 1.11 K. A modelfl^ that includes a small num- 
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ber of binding sites for ^He or atoms yields the func¬ 
tion form X = exp{a — R/T), where exp(a)/(l-|-exp(a)) is 
the minimum impurity concentration that blocks super¬ 
flux, and R includes the binding energy. This function 
form gives better fit (solid red line), with i? = 1.32 K and 
a = 2.18. The numbers reported in this paragraph are re¬ 
visions of those previously reported [l^ because they take 
into account the temperature scale revision mentioned 
earlier in this report. This energy value is higher than 
the predicted binding energy 0.7 K) of single 

^He atoms to dislocation cores. Although to our knowl¬ 
edge it has not been calculated, the binding energy to 
dislocation intersections should exceed this. These facts 
are consistent with the possibility that the flux extinction 
results from the ^He binding to dislocation intersections 
[S^l . where the ^He blocks the flux. To our knowledge, 
the binding energy for ^He at the intersections of dislo¬ 
cations has not been considered theoretically and is not 
known. 


It is perhaps useful to carry out quantitative estimates 
that relate to the decoration of such cores (or intersec¬ 
tions) with ^He atoms. There are a number of unknowns. 
One of these is the number of such structures that span 
the cell, Vycor to Vycor. We previously tookfl^ this 
number to be of order 10®. If we take the solid "'’He den¬ 
sity to be Ri 2.9 X 10^^ atoms/cm® we find that about 
6 X 10^ ^He atoms would be along a 2 cm direct strictly 
1-D pathway between the two Vycor rods; a pathway of 
diameter ~ 1 nm(^ would require ~ 5 x 10® atoms. Thus 
to fully decorate 10® such cores would require ~ 10^® 
atoms (or ^ 10®^ for a 1 nm diameter case). Since Cor- 
boz et al.fl^ have shown that the decoration of dislo¬ 
cation cores does not have to be complete to influence 
the superfluidity on the core, we rather arbitrarily take 
a smaller number, ^ 1 x 10®^ atoms. Of course, a much 
smaller number could also be relevant since only a local 
dense decoration would be needed along a short segment 
of a dislocation core or at an intersection to substan¬ 
tially reduce the flux. None the less, we continue with 
this number. Now, in our experimental geometry a ®He 
concentration of 10 ppm results in the presence of about 
5.4 X 10®^ ®He atoms in the cell, which, if uniformly dis¬ 
tributed is a ®He density of 2.9 x 10®^ atoms/cm®. 


So, we can ask from what volume the needed number of 
atoms (e.g. for the 10®® case) would have to diffuse in the 
solid to decorate the dislocation cores. If we assume that 
these atoms diffuse from a cylindrical region to a disloca¬ 
tion core on the axis of the cylinder we can estimate the 
radial distance from the core over which they would have 
to travel. We find this radial distance to be ~ 1.5 x 10“® 
cm. We can then ask how long this will take. Those 
from a region near the core will arrive relatively quickly 
while those from further away (an increasing number in 
any radial interval) will arrive later. For this we take the 
®He diffusion constant near 26 bar and 100 mK at ^ 10 
ppm (extrapolating from the work of Eselson et al. 3^) 
to be H ~ 1.5 x 10“® cm^/sec and < 6Dt, where 

t is the time, and we find the time to be about 2.5 x 10“® 


seconds. For higher ®He concentrations the atoms would 
have to travel less far, but for higher concentrations the 
diffusion constant is smaller. The time for flux recovery 
at Td is well documented in Fig. [SI but the time for flux 
extinction is not yet well determined; we can only say 
that it is apparently less than ~ 200 seconds. In spite 
of the approximations involved, if we compare the com¬ 
puted numbers for diffusion times with the times cited 
near Fig. [5] and Fig. [6] for the flux to be extinguished 
or recover, we conclude that diffusion to dislocations or 
their intersections is not likely a limiting factor in the 
flux change at T^; it would likely be a rather fast process 
for the concentrations that we have studied. 


3. A Helium Film 

Another unlikely possibility is the flow of a helium film 
along the surfaces of the Vycor and the experimental cell. 
Were this to be the cause of the mass flux, one might 
expect that the temperature dependence would behave 
like a Kosterlitz-Thouless transition. This behavior is not 
evident in the data shown in Fig. [Tjor other similar sets 
of data. None the less, were we to imagine a surface layer 
of liquid "'He at the walls with a thickness of two atomic 
layers (above a solid-like layer adjacent to most surfaces; 
above that layer "'He is known to be located adjacent to 
a wall in liquid mixture situations), the number of ^He 
atoms involved would be ^ 1.4 x 10®®. At 1 ppm ®He 
this is approximately the number of ®He atoms in the 
cell. Given diffusion times, it is conceivable such a film 
could be poisoned. But, the fact that annealing reduces 
or eliminates the ability of the solid to carry a flux when 
subsequently cooled argues strongly against a superfluid 
film as the carrier of the flux. 


4- Vycor Pore Openings 


Another possibility for what causes a reduction in the 
flux at Td is ®He accumulation at the openings to the 
Vycor pores. This possibility has recently been empha¬ 
sized by Cheng et al. 3^ . In their experiments a variation 
of our approach was used. Instead of a superfluid-solid- 
superfluid geometry, they used a solid-superfluid-solid ge¬ 
ometry. They were able to observe some temperature de¬ 
pendencies that are similar to those we have found in our 
various experiments, particularly the presence of Td- In 
their discussion of the pore opening scenario, which they 
supported by calculations of the temperature-dependent 
binding of ®He to various possible binding sites, the pic¬ 
ture is that at Td the ®He moves from the solid to the 
interfaces where the solid helium in the cell meets the 
liquid helium at the openings of the Vycor pores. 

For the flux to be fully blocked by ®He in this sce¬ 
nario, the openings of the pores where the superfluid in 
the Vycor meets the solid must be blocked by ®He. For 
our experimental apparatus each Vycor rod surface meets 
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the solid ^He over a macroscopic surface area of about 0.3 
cm^. Given the properties of Vycor, we estimate that the 
open area of this surface that is comprised of pore open¬ 
ings is no less than 0.084 cm^. Each pore of diameter 7 
nm will have an open area at the surface of the Vycor 
no less than 3.84 x 10“^^ cm^. The number of such pore 
openings is estimated to be at least 2.18 x 10^^. All of 
these need to be blocked by ^He. How much is needed at 
each pore opening? It is likely that 1 monolayer will not 
be adequate. We take as an estimate a distance along 
the pore of two pore diameters and presume that if this 
volume were to fill with ^He the pore would be blocked; a 
greater length certainly might be required. In the vicin¬ 
ity of 100 mK the expected phase separation for liquid 
^He indicates that if the ^He were to be in the pores and 
blocking the flux it would be due to a high concentra¬ 
tion normal mixture in the pores. Clearly our numbers 
provide only a rough estimate, but it is not unreasonable. 

To full such a volume for all of the pores would require 
Ri 2.6 X 10^® atoms. At 10 ppm, this is about 0.5 percent 
of all of the ^He available; for higher (lower) concentra¬ 
tions it is a proportionally smaller (larger) fraction. So, 
we can ask in this case about how long it will take for 
the ^He to accumulate at the pore openings. We take the 
same parameters of 26 bar, 10 ppm ^He with a diffusion 
constant of ~ 1.5 x 10“® cm^/sec and note that the ^He 
will have to travel macroscopic distances. For the case 
of 10 ppm ^He we find that the time required is ~ 20 
seconds. Changing the concentration changes the time 
required; increasing the length of the pore that needs to 
be filled with ^He increases it. These estimates are in 
very rough order of magnitude agreement with the times 
shown for recovery of the flux documented in Fig. [6l es¬ 
pecially when it is recognized that as ^He accumulates 
in the solid near a pore opening the local concentration 
may increase, which will cause the diffusion constant to 
decrease [13 . As we have noted, flux extinction is com¬ 
plete in ^ 200 seconds; recovery in ~ 500 seconds. 

To explore the flux increase for T < Td a.t the lowest 
concentrations it may be that the ^He is exhausted and 
incompletely effective in blocking the pores. For the case 
of 0.17 ppm ^He the filled cell will have 9.2 x 10^^ ^He 
atoms in it. As we have noted, to fill the pores to a depth 
of two pore diameters will require ~ 2.6x10^® atoms. We 
have seen that the flux recovers for ^He concentrations 
up to 10 ppm (in which case there are 5.4 x 10^^ ^He 
atoms in the cell). These estimates have a number of 
assumptions; in spite of them it is not clear why there is 
not enough ^He to completely block the flux. The tran¬ 
sition to low flux is prompt, but for low concentrations 
is not complete. One probable cause for this is that at 
a given temperature which is Td for the concentration in 
the solid, once ^He atoms begin to leave the solid ma¬ 
trix, the solid now contains a lower concentration of ^He. 
This naturally shifts the Td to a lower temperature and 
at the given temperature no additional ^He atoms leave 
the matrix. This observation explains the fact clearly 


seen in, for example. Fig. |4] and evident in other data 
sets, that the flux in the middle of the Td region is stable 
if the temperature is fixed. That is, the data taken while 
cooling or warming is stable once a fixed temperature is 
achieved. 


5. Behavior above and below Td 

It is not clear what might explain the universal tem¬ 
perature dependence above Td- And, unless the estimates 
made here are substantially in error, it is also not fully 
clear why the flux increases for low ^He concentrations 
below Td- One possibility below Td is that if the ^He 
has been exhausted, but the blockage incomplete, then 
a lower temperature could be expected to provide an in¬ 
crease in the superfluid density in the confining conduc¬ 
tance pathway. In very recent work we have reported ev¬ 
idence that for the region of this universal temperature 
dependence, T > Td, the limitation to the flux is unre¬ 
lated to the Vycor interface with the solid and takes place 
in the bulk solid [s^]. This suggests that this temperature 
dependence may be due to a temperature-dependent su¬ 
perfluid density along the conducting pathways. Addi¬ 
tional work is currently in progress and will be discussed 
more extensively in a future publication. 

IV. CONCLUSIONS 

We find the presence of ^He as an impurity in hep solid 
^He has a strong effect on the sharp flux reduction at a 
concentration-dependent characteristic temperature Td- 
On the other hand, we find that the presence of ^He does 
not alter the universal temperature dependence of the 
limiting mass flux above Td- The magnitude of the flux 
is typically sample-dependent, and sample-dependent at 
any concentration, but the temperature dependence is 
universal. The specific reason for the universal temper¬ 
ature dependence for T > Td remains unresolved. It is 
likely due to the physics associated with the conducting 
pathways. The results also suggest that the presence of 
^He does not destroy the apparent Luttinger-like behav¬ 
ior of the flux. More experimental and theoretical study 
of solid helium and ID superfluidity 0, H^l and its 
pressure dependence are needed. 
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